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^ ■ Abstract 



The stationary phase technique is used to calculate asymptotic 
formulae for S0(4) Relativistic Spin Networks. For the tetrahedral 



$—1 . 

spin network this gives the square of the Ponzano-Regge asymptotic 
formula for the SU(2) 6j-symbol. For the 4-simplex (lOj-symbol) the 
asymptotic formula is compared with numerical calculations of the 
Spin Network evaluation. Finally we discuss the asymptotics of the 



X 

; S0(3,l) lOj-symbol. 



1 Introduction 

Spin networks have been used to develop discrete models of quantum gravity 
called state sum models. The original state sum model due to Ponzano and 
Regge jlj used tetrahedral spin networks for the group SU(2) (6j-symbols) 
glued together to form a three-dimensional space-time manifold. The state 
sum model has a partition function which is determined by summing over 
the spin labels or parameters in the 6j-symbols, forming a discrete analogue 
of the functional integral for quantum gravity. These models were extended 
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to four-dimensional manifolds by considering spin networks based on a four- 
simplex ('lOj-symbol'), and simultaneously changing the group to SO (4), and 
later, the Lorentz group S0(3,l): the 'relativistic' spin networks [E]- 

A useful tool for analysing the physical content of these models is the 
analysis of the asymptotics of the 6j and lOj symbols. This involves deriving 
a formula for the behaviour of the values of the spin network when all the spin 
labels are simultaneously made large. Ponzano and Regge ^ conjectured a 
beautiful formula for the asymptotics of the 6j symbol expressed in terms of 
the geometry of the tetrahedron associated to the spin labels. This formula 
is essentially the Einstein-Hilbert action for the geometric tetrahedon, thus 
making an explicit connection between the state sum model and quantum 
gravity, at least in the asymptotic regime. The formula was later proved by 
Roberts |2] using the methods of geometric quantization. 

The asymptotic analysis for the lOj symbol was begun in [0] , which anal- 
ysed only some of the contributions in a stationary phase approximation 
(the non-degenerate ones, see below). Baez, Christensen and Egan ||8j jH] 
[Tn] performed some numerical calculations which gave a different scaling be- 
haviour, indicating that the remaining contributions (the degenerate ones) 
not analysed in [Hj must be important. 

In this paper we consider a fuller asymptotic analysis of the SO (4) lOj 
symbol, including the degenerate configurations, to give a complete picture of 
the asymptotics, thus explaining the numerical results in JUj. Some similar 
results to ours were recently obtained independently in [TT] . 

The paper develops some general machinery for evaluating the asymp- 
totics of relativistic spin networks. This machinery is applied to two specific 
cases: the tetrahedral network, and the 4-simplex (lOj symbol). The im- 
portance of considering the tetrahedal network is that the result is already 
known, since the relativistic tetrahedral spin network (for S0(4)) is just the 
square of the SU(2) 6j-symbol. Our asymptotic formula turns out to be ex- 
actly the square of the Ponzano-Regge asymptotic formula. Thus this gives 
essentially a simpler proof of the Ponzano-Regge formula than that given by 
Roberts. An interesting feature of our analysis, carried out using the station- 
ary phase approximation to an integral, is that it has both 'degenerate' and 
'non-degenerate' geometrical configurations. In this tetrahedral case, these 
two types of configurations have the same amplitude and decay at the same 
rate in the asymptotic region. 

Moving on to the SO (4) lOj symbol, we apply a similar analysis but 
find the degenerate configurations dominate in the asymptotic region with a 
slower decay rate; this explains the numerical results of [Hlj. We analyse the 
different types of possible stationary points in detail, giving the corresponding 
decay rates for the asymptotics. Some of the types of stationary point occur 
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for generic values of the spin labels on the 4-simplex while other types are 
non-generic. In some non- generic cases the degenerate configurations no 
longer dominate and we identify some of the other configurations which can 
dominate the asymptotic formula. 

The final section discusses the differences between the SO (4) (Euclidean) 
and the S0(3,l) (Lorcntzian) lOj-symbol. The stationary points of the Eu- 
clidean case all have analogues in the Lorentzian case; in particular there con- 
tinue to be non-degenerate stationary points for the Lorentzian lOj-symbol 
corresponding to geometric 4-simplexes in Minkowski space with spacelike 3- 
dimensional faces. However in the generic case the degenerate configurations 
still dominate the asymptotic formula. 
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2 Tetrahedron 

2.1 Evaluating SO(4) 6j-symbols 

We begin by calculating the evaluation of the Relativistic Spin Network 
Tetrahedron. The Relativistic Spin Network is the graph shown in figure ^ 
with edges labelled by half-integer spins. 
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Figure 1: Relativistic Spin Network Tetrahedron 

The evaluation of this network, a real number, is defined in terms of SU(2) 
spin networks by 




(1) 

In this equation the subscript R indicates a relativistic spin network, and 
SU(2) indicates an SU(2) spin network. Labels are omitted from the above 
equation but the edges on both SU(2) tetrahedra have the same labels as the 
relativistic one on the left hand side. The four theta symbols are formed by 
connecting each of the four vertices of the tetrahedron to themselves. 
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The definition of tlie SU(2) 6j-symbol gives 



{6j}s'(/(2) 




5(7(2) 
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SU{2) 
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SU(2) 




SU{2) 



(2) 



so tliat tlie SU(2) 6j-symbol is related to tlie relativistic spin network evalu- 
ation by 




SU{2) 



(3) 

Ponzano and Regge gave an asymptotic formula for the 6j-symbol ^J. 
Using notation {6j} for the value of the 6j-symbol and PR for the Ponzano- 
Regge approximation we write {6j} ~ PR to denote that {6j} is asymptot- 
ically equal to PR as all the spins are scaled upwards. Ponzano and Regge 
were not too specific about which scalings of the spins the formula should be 
valid for. For the formula to be useful it would be necessary to formulate this 
scaling in the most general possible way that is compatible with proving the 
result. However this paper is somewhat more preliminary as we are interested 
in establishing the asymptotic formulae in a restricted set of circumstances 
to give an outline of the asymptotic behaviour. The scaling that we use is 
to replace 2iki + 1 everywhere by a{2jki + 1) for fixed jki, allowing a — > oo 
through values where this makes sense (e.g. when a is an integer). 

The asymptotic formula means that lim^^oo({6j} - PR)a^l'^ = 0. The 
Ponzano-Regge formula is as follows 



PR 



V12ttV' 



cos 



E {^^^ + 1) 



Oki + -vr 



(4) 



where V is the volume of the tetrahedron with edge lengths J + |, and 9ki 
is the angle between the normal vectors to the faces k and /. 
Squaring this formula we obtain 



PR' 



^ (^l + cos (^{2j,i + 1)9^1 + 171^1 



(5) 



where V is the volume of the tetrahedron with edge lengths 2j + 1. 
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Notice that PR^ is composed of two parts, a cosine term and a constant 
term, both of which have the same amphtude. We call the cosine term 
[PB?)cos and the constant term {PR^)const so that 

PR' = {PR'),onst + {PR')cos. (6) 

The evaluation of the relativistic spin network can also be computed using 
the integral formula |1] 

/= (_l)E.<,2iu /■ T\Tr{pMixkXi')dx (7) 

JxeSU{2r k<l 

where pki is the representation of SU(2) on edge kl and Xk is an element of 
SU(2). We integrate over four x's, one at each vertex of the tetrahedral spin 
network. The integration uses the Haar measure, normalised to unity. We 
calculate /' = I{-l)^k<i'^jki go that /' = {Qj^. 

In the next sections we will rewrite the integral in such a way that the 
stationary phase formula can be used. Each stationary phase point is a 
certain configuration for a geometrical tetrahedron (or a generalisation of 
this). We will show how the different configurations give the two terms in 
equation 



2.2 Kirillov Character Formula 

We begin by rewriting ((7j) using the Kirillov character formula. 



sin(rfc;) 4,,es2 47r 



where ^ki is an element of the Lie algebra of SU(2) defined by 

exp(^M) = xixl^, (9) 

and Tki = \Cki\- The ambiguity in the definition of ^ is fixed by requiring 
^ki = \(,ki\ to be the angle between Xk and xi thought of as vectors in = 
SU(2). This gives a unique ^ for all angles < r < vr. The case r = vr, when 
two x's are anti-parallel, requires special treatment, as the Kirillov formula 
does not work directly there. One could use a simple modification of the 
Kirillov formula which would work in a neighbourhood of r = vr. However 
consideration of the cases in which r = vr can be effectively bypassed, as 
discussed below. 
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Carrying out the integral in equation (jH)) gives the Weyl character formula 
used by Barrett and Williams to find the non-degenerate stationary points. 

rrv/ / -ixN sin((2jfc; + l)rki) 

Trip,,{x,x, )) = -^^ (10) 

It is important here to use the Kirillov formula, rather than the Weyl formula, 
as it enables us to evaluate the contribution from geometries where Vki = 0. 




Figure 2: Interpretation of the Kirillov Character formula 

Considering elements Xk of SU(2) as vectors in S^, and replacing 2jM + 1 
by a{2jki + 1) to examine the asymptotics when a — oo, we obtain 



(47r)6(27r2)4 ^ ^ 

r r f -|-|- a{2jki + l)rki \ ^iay:,.^A2ju+i)iu..vM^y^^ 

7,6(53)4 4g(s2)« sin(r«) ) 

We call the exponent J2k<ii'^jki + ^)iki-yki in equation (|TT|l the action, and 
denote it by S. 

We consider the integration to be performed in two parts, firstly fixing 
Xi to be at the north pole of S*^, X2 to be in a plane containing xi, and X3 
to be in a 3 dimensional hyperplane containing the plane. The second part 
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is to integrate over the symmetry of the integrand (which only depends on 
relative angles), rotating all the x's and y^s at once. 

The symmetry group SO (4) is 6 dimensional, the original integration was 
over 24 dimensions, so the first part of the integral is over 18 dimensions. 
Integrating over the symmetries will multiply the result by a volume factor. 

2.3 Method of Stationary Phase 

We use the method of stationary phase to calculate the asymptotic contri- 
butions to the first part of the integral. The general stationary phase formula 
is 

(12) 

where H{x) is the Hessian matrix for 0(x), and sgn(if) is the number of 
positive eigenvalues minus the number of negative eigenvalues of H. 

Writing for the asymptotic approximation to /' i.e. I' ~ J^, this gives 

(13) 

where 2^7r^ is the volume of the space we have quotiented out of the integral 
to remove the SO (4) symmetry. The assumption is that in the remaining 
integral the stationary points are discrete. Notice that all terms in this 
formula will scale as a~^. 

We must find the stationary points of the action 

S = Y,{2jui + mum (14) 

k<l 

with respect to varying Xj and yku remembering that each S,ki depends on the 
x's via equation ©. 

There are a number of types of solutions for which the action is stationary. 
Firstly the solutions can be classified by the dimension of the linear subspace 
which is spanned by the unit vectors Xj. Secondly, there is a qualitative 
difference between solutions for which no pair of these vectors are parallel, 
i.e. Xi 7^ ±Xj for all i ^ j-, and solutions in which at least one pair are 
parallel. The former will be called 'non-parallel' solutions and the latter 
'partly parallel' if not all Xi are parallel, and 'totally parallel' if Xi = ±Xj for 
all i,j. Both of these classifications measure the amount of degeneracy of 
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the solutions, and there is some relation between them. In the two examples 
considered below, the highest possible dimension solutions for a given spin 
network are necessarily non-parallel. These solutions will be called 'non- 
degenerate'. At the other end of the scale, the totally parallel solutions are 
clearly the same as the one- dimensional solutions. 

For the parallel solutions, one only needs consider the cases when Xi = xj. 
Replacing Xk by —Xk has a very simple effect on the value of the integrand 
in (jZj): it is multiphed by the factor 

The factor is equal to 1 if the parity admissibility condition is satisified, and 
— 1 otherwise. This means that if the admissibility condition is not satisfied, 
then the contributions from stationary points at Xk and —Xk always cancel 
and so the asymptotic formula always gives zero. If the condition is satisfied 
then one only need consider one of these two possibilities, and one can always 
choose Xi = Xj rather than Xi = —xj. Thus one never needs to do a detailed 
calculation of the contribution of stationary points for which the angle r = vr. 



2.4 Results 

The following sections show that whenever the spins are the edge lengths of 
a Euclidean tetrahedron the stationary phase formula has discrete stationary 
points which are one of two types, either non-degenerate, or one-dimensional. 
Accordingly, the stationary phase formula can be written 

The contribution from non-degenerate configurations, Cnon-deg, gives exactly 
the oscillatory part of the Ponzano-Regge squared formula 

The contribution from one-dimensional configurations, Ci_rf, gives exactly 
the constant part of the Ponzano-Regge squared formula 

= iPR^)const. (16) 

Putting these two results together, this shows that the asymptotic formula 
for the relativistic spin network tetrahedron agrees exactly with the square 
of the Ponzano-Regge formula whenever we have a tetrahedron which can be 
embedded into 3-dimensional Euclidean space, 

/: = PR'. (17) 
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2.5 Non-parallel configurations 

Varying the action with respect to the ?/'s we find that the action is stationary 
when i/ki is parallel or antiparallel to C,ki- In other words, 

Vki = efc/]|^, (18) 

iQkll 

with eki = ±1. 

Varying the action with respect to the x's, and using a Lagrange multiplier 
Aj for the constraint Xj.Xj = 1 we obtain four equations, one for each face of 
a tetrahedron. 

ei2(2jl2 + 1)X2 , ei3(2ji3 + 1)X3 ei4(2ji4 + 1)X4 _ 

r-7 r \ ^ \ r-7 ^ — AlZi 

sin(ri2) sin(ri3) sin(ri4) 

ei2(2jl2 + ^ e23(2j23 + 1)3:3 ^ e24(2j24 + l)a^4 _ 

sin(ri2) sin(r23) sin(r24) ^ ^ , . 

ei3(2ji3 + 1)3^1 , e23(2j23 + 1)3^2 , 634(2^34 + l)a;4 _ , 

r-7 ^ \ r-7 ^ \ r-7 ^ — A3a;3 

sin(ri3) sin(r23) sin(r34) 

ei4(2ji4 + l)xi ^ e24(2j24 + 1)3^2 ^ e34(2j34 + l)a^3 _ ^^^^ 
sin(ri4) sin(r24) sin(r34) 

The equations show that there is at least one linear relation between the 
x's and so require the x's to lie in a 3- dimensional hyperplane. Taking v^i to 
be the unit vector in the direction Xk x xi (using the vector cross product in 
the 3d hyperplane) gives 



ei2(2ji2 + l)vi2 + ei3(2ji3 + 1)^^13 + ei4(2ji4 + l)vu = 

-ei2(2jl2 + 1)^12 + e23(2j23 + l)t'23 + e24(2j24 + 1)^^24 = 

-ei3(2jl3 + 1)^^13 - e23(2j23 + 1)^^23 + e34(2j34 + 1)^34 = 

-ei4(2ji4 + 1)V14 - e24(2j24 + 1)^24 " e34(2j34 + 1)^34 = 



(20) 



2.5.1 Three-dimensional configurations 

Now consider the case of a three-dimensional solution. Setting Vki = eki{2jki+ 
l)vki, the above equations show that the Vki are the edge vectors of a tetra- 
hedron with edge lengths 2jki + 1. Each equation describes the face of a 
geometrical tetrahedron as shown in figure El Notice that the four equations 
are not linearly independent - the sum of the first three gives the fourth 
equation. 

Each normal vector to a face of the tetrahedron; it could be inward 

or outward pointing. Let be the outward normal to the i-th face. Then 
Xi = aiUi with Oj = ±L 
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Figure 3: Non-degenerate stationary point 



In figure El the edge vectors for tlie tetrafiedron are given by tlie expression 
^{'^jki + l)'^fc ^ ni = ±(2jfci + VjakdiVkh witli ± indicating one overall choice 
of sign. This means that, for figure El the e are determined by eui = Xauo-u 
with A = ±1. This means that many choices of eki give no solutions at all. 

Every stationary point is of this form. Consider any set of x's and e's 
satisfying equations ((201), then this corresponds to a geometric tetrahedron. 
By acting on the x vectors with an element of SO (4) this tetrahedron can be 
rotated to correspond to the tetrahedron drawn in figure El possibly with all 
of the arrows indicating the directions of the V reversed. 

Given one stationary point {xj, eki\ we can generate all of the others using 
combinations of two operations. The first operation is to swap the sign of 
one or more of the x's. Since Xi is fixed at the north pole we are not at 
liberty to swap its sign. For = ±1, 



Xi Xi 

X2 a2X2 

X3 0-3X3 

X4 — > (T4X4 

Cfc/ —>■ o'kCTieki- 



The second operation is to simultaneously swap the sign of all of the e 
whilst leaving the x unchanged 
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Xi ^ Xi 
^kl — ^ —^kl 

'^vertex ^ki IS an integer, each of the 8 stationary points related by the 
first operation gives the same contribution to the integraL The stationary 
points which are related to these by the second operation give the complex 
conjugate contribution. 

If Tlivertexiki is not an integer (i.e. is an odd half-integer), the stationary 
points for the x's cancel and there is no contribution from the non-degenerate 
stationary points. This corresponds to the admissibility condition for spin 
networks. For a spin network to be called admissible the j's labelling the 
edges incident to each vertex must sum to an integer (and satisfy triangle 
inequahties) . The admissible spin networks are the only ones with non- 
zero evaluation. If the admissibility condition is not met there will be no 
stationary point contribution. The stationary phase calculations also give the 
triangle inequalities for spin networks, since to be able to form the triangular 
faces of a tetrahedron each spin in a face must not be greater than the sum 
of the other two. 

The Hessian is an 18 x 18 matrix with components Hki — q'^J).^^ where 
Ui - . . Mis are the coordinates we are integrating over. For any particular 6j- 
symbol the eigenvalues of the Hessian can be found using computer algebra. 
The Hessian in the e's all positive case has 2 positive eigenvalues and 16 
negative eigenvalues. Hence sgn(if) = —14. The e's all negative case has 
sgn{H) = 14. 

These contributions will combine to give an asymptotic expression of the 
form 

= ^ cos ( Y,(2jki + + 7r/2 ] (21) 

\k<l / 

The determinant of the Hessian can be evaluated using computer algebra 
for any particular 6j-symbol. The amplitude of the contribution has been 
calculated for a large number of cases and compared with the amplitude of 
the cosine term in the square of the Ponzano Regge formula. Examining a 
large number of cases, we find that they are always equal: 
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More precisely, the amplitudes agree exactly for all Euclidean tetraliedra, 
which are refered to as type I in Ponzano-Regge pp. Ponzano and Regge's 
type II and III tetrahedra, the Lorentzian and transition cases do not have 
an analogue here. However we have been unable to prove this formula in 
general, as the algebraic expressions for the determinant of the Hessian are 
too large. 

The final result is 

Cnon-deg = COS | ^(2jH + l)rkl + 7r/2 j . (23) 

2.5.2 Non-parallel 2-dimensional configurations 

The only lower dimensional non-parallel configurations possible are where all 
the x's lie in a 2 dimensional hyperplane. If the admissibility conditions for 
the tetrahedral spin network are satisfied then there are no 2 dimensional 
solutions to the stationary phase equations (j20|) . 

If the admissibility conditions are not satisfied then the contributions for 
different e's exactly cancel. Hence the lower dimensional configurations do 
not contribute to the asymptotic formula. 



2.6 One-dimensional configurations 

The one- dimensional configurations involves all x's parallel or antiparallel. 
As noted above, one need only consider the case when all x's are actually 
equal, the other cases giving the same asymptotics. 

Returning to the action ()14|) if we vary the y's the action is stationary in 
all cases. If we vary the x's we find that for the action to be stationary 

(2ji2 + i)yi2 + (2ji3 + i)z/i3 + (2ji4 + i)yu 

-(2ji2 + 1)2/12 + (2j23 + 1)2/23 + (2j24 + 1)2/24 
-(2ji3 + l)yi3 - (2j23 + 1)2/23 + (2j34 + 1)2/34 
-(2jl4 + 1)2/14 - (2j24 + 1)2/24 - (2j34 + 1)2/34 

This means the 2/'s are parallel to the edge vectors of the same tetrahedron 
we found in the non-degenerate case. Indeed, setting Yki = {2jM + l)yki we 
now have the Y^i as edge vectors of a tetrahedron with edge lengths 2j -|- 1, 
as in figure |3 (or with the direction of the arrows swapped) 

There will only be a solution to the stationary phase equations if the j's 
satisfy the triangle inequalities. There are 2^ stationary points of this type 
for the x's (since we can choose each of X2 . . . X4 to be parallel or antiparallel 










(24) 
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Figure 4: Degenerate stationary point 



to Xi) and 2 sets of y's for each (the vectors could be as drawn in figure 
0] or reversed). The contributions are identical if the parity admissibility 
condition is satisfied, and exactly cancel if not. 

The Hessian for the all x's parallel case has 9 positive and 9 negative 
eigenvalues. Hence sgn(if) = 0. These contributions will be non-oscillating, 
since S = for one-dimensional configurations. 

The amplitude of the contribution has been calculated for a number of 
cases, and compared with the amplitude of the constant term in the square 
of the Ponzano Regge formula. Again we find that in all cases 



Denoting the contribution from one- dimensional configurations by C 



l-d 



Ci-d = (26) 

2.7 Higher-dimensional parallel configurations 

There is never any contribution from partly parallel configurations where 
some but not all of the x's are parallel. Again if the admissibility conditions 
for the tetrahedral spin network are satisfied then there are no solutions 
to the stationary phase equations. If the admissibility conditions are not 
satisfied then the contributions for different e's exactly cancel. 
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3 4-Simplex 

3.1 Evaluating SO (4) lOj-symbols 

We now repeat the previous calculation for the Relativistic Spin Network 
4-simplex, or SO (4) lOj-symbol. 




Figure 5: Relativistic Spin Network 4 Simplex 



We use the same integral formula, now with 5 vertices and 10 edges. 

I ^ (_i)E.<,2,., f T\Tr{pM{xkX^')dx (27) 

Set /' = (— l)Sfc<i2.7*!j simplify the equation. Replacing 2jki + 1 by 
a(2jki + 1) to investigate the asymptotics, and using the Kirillov character 
formula we obtain 



^ (47r)io(27r2)5 ^^^^ 

For the lOj symbol many of the stationary points are no longer points 
but higher dimensional manifolds. We integrate over this space of stationary 
points as for the symmetry group in chapter 2. Each degree of freedom (i.e. 
dimension) in the stationary phase solution manifold reduces the size of the 
Hessian by 1, and hence changes the scaling behaviour of the asymptotic 
contribution by a 2. We are integrating over 35 dimensions (15 for the x 
variables and 20 for the y variables), and there are always 6 symmetry degrees 
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of freedom from 5*0 (4). If there are an additional n degrees of freedom in 
the stationary phase solution then the amplitude of the contribution will be 

35 — 6 — Ti -1 fi 9 — n . , 

a —a^^ = a-— (29) 

where the factor of a^^ comes from the product in the integral. 
The stationary points are the stationary points of the action 

S = Y,{2jki + mu ■ Vki. (30) 

k<l 

We consider the contribution from non-degenerate, lower dimensional, 
and parallel configurations. 



3.2 Results 

The following sections show that the stationary points in the asymptotic for- 
mula are one of a number of types. For a generic set of spins the scaling be- 
haviour of these types follows. The one-dimensional stationary points 
give a contribution which scales as a"^. The four- and three-dimensional 
non-parallel stationary points give contributions which scale as al- 
though it is possible that in some cases one or both of these types may be 
absent (due to inequalities which the spins must satisfy). The partly par- 
allel stationary points have only one generic contribution (again subject to 
inequalities). This is a three-dimensional stationary point with just two of 
the X variables parallel. Again this scales as 

This means that that asymptotically, the evaluation of the lOj symbol 
is dominated by the degenerate one- dimensional stationary points, whose 
contribution in general decays constant times a , 

I'a ~ Koi-\ (31) 

This is in contrast to the 6j symbol case where the contribution from the 
degenerate and non-degenerate stationary points had equal amplitude. 

The are some non-generic cases of spins where this behaviour can differ, 
as mentioned in the text below. For example, at the limit of admissibility, 
the scaling of the one-dimensional configurations can increase to a"^/^ or 
this type may be absent; in these cases it is expected that the dominant 
contribution will come from a non-generic partly-parallel configuration. For 
one or two vertices at the limit of admissibility, these configurations give an 
asymptotic contribution which scales as . This agrees with the numerical 
results obtained by Baez, Christensen and Egan[TU]. Other examples of non- 
generic behaviour are detailed in the following sections. 
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3.3 Non-degenerate 4 dimensional configurations 

First we consider non-degenerate solutions where the Xi have a 4 dimensional 
span, following W\. 

Varying the action with respect to the y's we find that the action is 
stationary when y^i is parallel or antiparallel to ^ki- Let e^/ = ±1 then 

Vki = (^kiTT-^ (32) 
The action at the stationary points in the y's becomes 

S = Y,^u{23M + l)rki (33) 

k<l 

A set of 5 non-degenerate vectors with a 4 dimensional span determine a geo- 
metric 4-simplex up to scale P|. For a geometric 4-simplex we have Schlafii's 
identity 

^Afc,rf0« = O (34) 

k<l 

where are the angles between the outward normal vectors nk-,ni to the 
tetrahedra k and / of the 4 simplex. The A^i are the areas of the triangles of 
a geometric 4-simplex determined by the (pki (up to overall scaling of the 4 
simplex) . Varying the action, using a Lagrange multiplier A for the constraint 
equation (jSl)), gives 

dS = ^ eki{2jki + 1) drki = A ^ dcpki (35) 

k<l k<l 

If we take the x's to be the outward normals to the tetrahedra of the 
4 simplex then Vki = 4>ki- If an x is swapped to be an inward normal then 
f^ki = — <Pki'i hence dvki = —d(j)ki for all r's involving the inward x. 

Taking to be the outward normal, and introducing variables which 
are +1 if is outward and —1 if is inward pointing, we have Xk = dk^k- 
For each triangle we obtain 

eH(2jH + 1) = XataiAki (36) 

Taking A = ±1 to fix the overall scale of the A's we see that A^i = (2jfc;+l) 
and eki = Xa^ai for each A;, /. Hence at a non-degenerate stationary point the 
Tki are the angles between the inward or outward pointing normals to a 4 
simplex with areas (2j + 1). 

If the parity admissibility conditions (sum of j''s at every vertex is an 
integer) are satisfied for the lOj-symbol, each solution for A = +1 gives the 
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same contribution, and each A = — 1 solution gives the complex conjugate 
contribution. These terms will combine to give a cosine contribution to the 
integral. There will be 2^ distinct stationary points for the x's (y because 
there are 5 a's which can each be ±1 but swapping the sign of them all 
corresponds to — / G S0(4)) and 2 possible sets of e's for each, coming from 
the choice of A. Changing A changes the sign of all of the e's, swapping the 
direction of all of the yu- 

If the parity admissibility conditions are not satisfied the terms will cancel 
giving no net contribution. Inequalities analogous to the triangle inequalities 
for vertices in the 6j-symbol will occur in this case because to be able to 
form a tetrahedron with four triangles of area (2j + 1), each area must not 
be greater than the sum of the other three. 

There may in general be more that one geometric 4 simplex with a given 
set of areas, for example see Tuckey's example in which is a Euclidean 
4-simplex for t < 8/3. When this occurs each possible configuration will 
have a contribution to the asymptotic formula. From one solution to the 
stationary phase equations swapping and will yield another solution 
[TTj . In general these are the only pair of vectors which can be interchanged 
like this since xi . . .x^ are constrained to lie in hyperplanes. 

Terms from the Hessian and the integrand show that the amplitude of 

9 

the contribution from the non-degenerate stationary points will scale as a~^. 
3.3.1 Example: Regular 4-simplex 

For the regular 4 simplex with a(2jfc; + 1) = f3 for all k, I we get a contribution 
Cnon-deg from nou-dcgenerate geometric 4-simplexes. A corollary of Bang's 
theorem |7j tells us that there is only one geometric 4 simplex with all areas 
equal . The corresponding contribution has been calculated by computing 
the value of the non-exponential part of the integrand and the Hessian at the 
stationary points. 

C„„„_., = M^T^ ™^ (lO/3cos-(-l) + (37) 
3.4 One-dimensional configurations 

For this case we take all the x's to be parallel or antiparallel. Each of the 
2^ cases has an identical contribution to the integral. The action is always 
stationary with respect to varying the y's in this case. Varying the action 
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with respect to the x's gives five equations in the y variables. 



(2ji2 + l)yi2 + (2ji3 + I)l/13 + (2jl4 + I)l/14 + (2jl5 + l)l/l5 = 
-(2ji2 + l)yi2 + (2j23 + l)l/23 + (2j24 + 1)^/24 + (2j25 + l)Z/25 = 
-(2ji3 + l)yi3 - (2j23 + 1)2/23 + (2j34 + 1)2/34 + (2j35 + 1)2/35 = 
-(2ji4 + 1)2/14 - (2j24 + 1)2/24 - (2j34 + 1)2/34 + (2j45 + 1)2/45 = 
-(2jl5 + 1)2/15 - (2j25 + 1)2/25 - (2j35 + 1)2/35 " (2j45 + 1)2/45 = 






(3^ 







Setting Yki = {2jM + l)yki these equations correspond to the vector dia- 
gram in figureini The front, top, bottom, left and right faces each correspond 
to a stationary phase equation above. The back face is added in to complete 
the figure. 




^25 



Figure 6: 4-simplex totally parallel stationary phase point 

Generally the y vectors subject to the above set of equations will have 
8 remaining degrees of freedom. This is because there are 20 y variables 
and 12 independent equations. These 8 degrees of freedom consist of 3 from 
rotations in S0(3) and 5 deformations of the above 'cube'. The angles in the 
cube in figure IHl can be varied as long as the edge lengths remain unchanged. 

To see that there are 5 degrees of freedom in deforming the 'cube' notice 
that there is one angle to determine the geometry of the back parallelogram 
face, then 2 spherical coordinate angles each will determine the directions of 
two diagonally opposite vectors between the back face and the front face (2/14 
and 2/23 say). Once these choices have been made the rest of the vectors are 
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fixed up to a finite set of discrete possibilities, by requiring the edge lengths 
to be 2jki + 1. There will be certain inequalities which must be satisfied to 
ensure that the figure can be embedded in Euclidean space. Thus the parallel 
stationary point contribution is of the form a~~ = 

For a stationary phase solution of this type to exist and for the contribu- 
tions from parallel and antiparallel terms to add up requires satisfaction of 
the admissibility conditions for the lOj-symbol. 

These results agree with numerical calculations of lOj symbols carried 
out by Baez, Christensen and Egan. In and ^U] numerical values of the 
lOj symbols are calculated and discussed. They see lOj symbols decaying as 
0(«-2). 



3.5 Partly parallel configurations 

We must consider the case where some but not all of the x's are parallel. 
There are a number of ways this may happen as there can be one or two 
directions each parallel to two or more vectors. Unlike the 6j-symbol, for 
the lOj-symbol contributions may come from configurations of these types, 
depending on the labelling by spins. 

The only generic case is when just two of the x vectors are parallel. Then 
the X span 3 dimensions and are the normal vectors to a tetrahedron. For 
this stationary point to exist inequalities between the spins must be satis- 
fied so that the tetrahedron may be embedded in 3 dimensional Euclidean 
space. The x and y vectors are uniquely determined by the stationary phase 
equations so there are no degrees of freedom. The contribution from this 

9 

configuration will scale as a~^. 

Non-generically, with two vectors parallel, if the spins satisfy a set of 
particular equalities, it is possible for one of the stationary phase equations 
to become trivial, allowing up to 4 degrees of freedom. For xi and X2 parallel 
this non-generic stationary point occurs only if 



Jl3 
jl5 



J23 
hi 



J45- 



(39) 



This gives a contribution scaling as a~^. 

If two pairs of x vectors are separately parallel there can be a stationary 
point only if the spins satisfy certain equalities. If Xi and X2 are parallel, 
and X3 and x^ are separately parallel then this stationary point occurs only 
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if there exists a set of e such that the j's satisfy 

335 = J45 (40) 

ei3(2ii3 + 1) + ei4(2ii4 + 1) + €23(2^23 + 1) + £24(2^24 + 1) = 0. 

The X vectors span 3 dimensions and have up to 3 degrees of freedom, whilst 
the y's are fixed. The contribution from this configuration, when it occurs, 
will scale as a~^. 

Three x vectors parallel can be a stationary point only as a special case 
of the two-dimensional configuration. The spins must satisfy a number of 
equalities so this is not a generic contribution. If there is no 2 dimensional 
stationary point with the x's non-parallel then the contribution from this 
configuaration when it exists will be of order a" 2. If there is a 2 dimensional 
non-parallel stationary point then this configuration will just be part of that 
stationary point. The |/'s are uniquely determined by the stationary phase 
equations. The same happens if three x vectors are parallel and the remaining 
pair of X vectors are separately parallel, however the contribution when there 
is no 2 dimensional stationary point will be of order a"^. 

Finally, four of the x vectors may be parallel. This is a non-generic 
stationary point since there is an equality between the spins. If Xi,X2, x^, X4 
are parallel then this stationary point occurs only if there exist a set of e such 
that 

ei5(2il5 + 1) + e25(2j25 + 1) + e35(2j35 + 1) + e45(2j45 + 1) = 0. (41) 

There are up to 2 degrees of freedom in the y variables and one degree of 
freedom in the x variables, so the contribution from this stationary point will 
generally scale as a'^. 

3.6 Non-parallel lower dimensional configurations 

There may also be stationary phase points for configurations in which the Xi 
are not parallel but lie in a hyperplane of dimension less than 4. 

3.6.1 3 dimensional configurations 

Consider configurations where Xi lie in a 3 dimensional hyperplane but no 
pairs are parallel. 

Using the vector cross product in the 3 dimensional hyperplane we write 
Vki for the unit vector in the direction x^ x xi, and V^i — eki{2jki + i)vki- The 
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stationary phase equations show that yki = tu^ki and 



-Vi2 + V'23 + V2A + V^25 
-^13 - V^23 + V^34 + ^35 
— Vi4 — V24 — V34 + V45 

-^15 - V^25 - V^35 - ^45 























(42) 



Notice that these equations are very similar to those obtained for the 
totally parallel case in ()38p. However since Vki is orthogonal to Xk for each 
/ we see that each equation in ()42|1 relates vectors which lie in a plane. The 
geometric figure differs from the totally parallel case because here all its faces 
are planar. The geometry is more constrained than in the totally parallel 
case. 

The geometric interpretation of a solution to these equations is a degen- 
erate case of a 4 simplex in which every triangle shares a common direction. 

The figure in the totally parallel case had 5 degrees of freedom. Now we 
are looking at the same figure with the additional constraint that each of 5 
faces must be planar. Therefore we expect that when solutions exist to these 

equations there will be no further degrees of freedom available, and that 

9 

the contribution to the integral will be a~2 multiplied by some oscillating 
function of modulus 1. 

3.6.2 2 dimensional configurations 

Consider configurations where not parallel but lie in a 2 dimensional 

hyperplane. The unit bivectors formed from the x-i will be equal (up to sign). 
Write 5ki for the sign of Xk A xi relative to xi A X2 (5i2 = !)• 

Then the stationary phase equations show that Uki = ^ki^^\ and 



ei2(2ji2 + l)5i2 + ei3(2ji3 + 1)^13 + ei4(2ji4 + l)5i4 + ei5(2ji5 + 1)^15 = 

-ei2(2jl2 + 1)^12 + e23(2j23 + 1)^23 + £24(2^24 + 1)^24 + e25(2j25 + 1)^25 = 

-ei3(2jl3 + l)'5l3 - e23(2j23 + 1)'^23 + e34(2j34 + 1)'^34 + e35(2j35 + 1)'^35 = 

-ei4(2ji4 + l)5i4 - e24(2j24 + 1)^24 - e34(2j34 + 1)^34 + e45(2j45 + 1)^45 = 

-ei5(2jl5 + 1)^15 - e25(2j25 + 1)^25 " e35(2j35 + 1)^35 " e45(2j45 + 1)^45 = 



These configurations will contribute to the integral whenever the above 
set of equations can be satified. This only happens for special values of j; 
generic lOj-symbols will not have stationary points of this type. 



(43) 
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For example of when this type of stationary point does occur consider 
the regular lOj-symbol with all j's equal. Arrange the x vectors in a plane so 
that moving clockwise round the vectors they are in ascending order. Then 
all bi vectors are equal and every 6 is +1. 

Choosing eis = —1, = —1, 625 = —1, £35 = —1, and all other e's equal 
to +1 gives one solution to the stationary phase equations. 

We have four degrees of freedom to choose the angles from xi to the 
other four Xi in the plane, and for all of these configurations the action will 
be stationary. The action will be zero for all of these configurations, so the 
contribution to the integral will be a constant multiplied by a~2. 

3.7 Limit of admissibility 

A vertex of a spin network is said to be at the limit of admissibility if the 
spins labelling the edges incident to it only just satisfy one of the admissibility 
inequalities for that vertex. For a vertex of a 4-simplex one of the spins is 
equal to the sum of the other three, for example 

jl2 = jl3 + jl4 + jl5- (44) 

At the limit of admissibility the behaviour of the stationary point equa- 
tions is typically different to the generic behaviour outlined above. However 
with the scaling given by 2j + 1 = a{2J + 1) for some constant J, this is 
not apparent because the scaling does not preserve the condition (jH)). Thus 
the only way to investigate the behaviour at the limit of admissibility in the 
asymptotic limit is to use a different scaling, namely 

j = aJ 

with constant J and a — > 00. The effect this has on the stationary phase 
formalism is to change the action (the part multiplying a) to 

S = Y,i:^3ki)iki-yku (45) 

k<l 

instead of (jHHj). Consequently the stationary point equations are the same 
but with (2j + 1) everywhere replaced with 2j. The following considera- 
tions indicate how the above types of stationary point behave at the limit of 
admissibility with this alternative scaling. 

The asymptotics of the one- dimensional configurations is modified in this 
situation. If one or more vertices of the lOj symbol are at the limit of ad- 
missibility the equation corresponding to that vertex will require the y's to 
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be parallel. If only one vertex is at the limit of admissibility (see figure 
El) then the degrees of freedom available in deforming the cube will drop 
in the asymptotic limit from 5 to 2. The contribution in this case will be 

9 — 2 7 

a 2~ = a" 2, and may no longer be the leading contribution to the limit. 




Figure 7: 4-simplex totally parallel stationary phase point - vertex 1 at limit 
of admissibility 

If more than one vertex is at the limit of admissibility there will in general 
be no solutions to the stationary phase equations with all x parallel. 

For the partly-parallel configuration with x\ and X2 parallel the equations 
()39|) show that if vertex 1 is at the limit of admissibility then vertex 2 must 
also be. This stationary point is expected to give the leading-order contri- 
bution to the asymptotic formula of in the case where two vertices are 
at the limit of admissibility. This agrees with the numerical results in JT] . 

If two pairs of x vectors are separately parallel, the equations are com- 
patible with two vertices at the limit of admissibility but the contribution 

7 

decays as and so this case does not give the leading contribution to the 
asymptotic formula. 

With four X vectors parallel, the stationary equations (j^T|) become 

ei5jl5 + £25^25 + £35^35 + £45^45 = (46) 

with this scaling. Some choices of e reproduce the limit of admissibility 
condition for vertex 5. This gives a term of order in the case of one 
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vertex at the limit of admissibility. This is the leading order contribution, 
again agreeing with the numerical evidence in ^T] . The contribution will die 
off more quickly if more than one vertex is at the limit of admissibility. 
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4 Lorentzian lOj Symbol 

4.1 Evaluating SO(3,l) lOj-symbols 

In this section the relativistic spin networks for the Lorentz group S0(3, 1) 
are analysed in a similar fashion to the Euclidean case. In the definition 
one replaces constructions related to SO (4) with the analogous concepts for 
S0(3, 1). In particular, the spin labels for the Lorentz group are non-negative 
real numbers. The label for the kl-th edge is denoted pki- There is an inte- 
gral formula for S0(3, 1) spin networks called the regularized evaluation, 
which takes the following form for the lOj-symbol (5 vertices and 10 edges). 
This integration is carried out over Hyperbolic space (the three dimen- 
sional subspace x.x = — 1 of — h ++ signature Minkowski space), instead 
of 5*^ for the SO (4) case studied previously. The distance rki between two 
variables Xk,xi G is defined by coshr^; = —Xk ■ xi. The evaluation is 



Again there is a vector x for each vertex of the spin network, and X2 ■ ■ - x^ are 
integration variables. For this formula, Xi is fixed at (1,0,0,0), removing a 
symmetry which would otherwise multiply the integral by the infinite volume 
of hyperbolic space. Note that the formula in ^2] for the evaluation has 
been multiplied by Y[k<iPki so that the analogy with the previous sections 
is clearer. Given two points Xk and xi in H^, it is possible to define a vector 
tangent to Xk which indicates the position of xi. Define to be the initial 
velocity of the geodesic •yit) which has endpoints 7(0) = Xk and 7(1) = xi. 
Then \^ki\ = ^ku the distance between Xk and xi. With this definition, 
^ki £ ^-f^zfe- However it is possible to identify all of these tangent spaces if 
desired using an appropriate trivialisation of the tangent bundle. 
Using the Kirillov formula 




(47) 




and replacing pki by apki to effect the scaling gives 



1 



(27r2)4(47r) 
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k<l 



akiPkirki 
sinh(rfc/) 
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We use the stationary phase method as before to examine the asymptotics 
as q; ^ oo. It is necessary to find the stationary points of the action S — 
YlikKiPkii^ki-yki) with respect to variations in x and y. 

4.2 Results 

The stationary points can be classified in the same way as for the SO (4) 
case. The stationary phase equations for S0(3, 1) are more or less exactly 
the same as for SO (4), except the vectors are in Minkowski space rather than 
four-dimensional Euclidean space. One qualitative difference is that the Xi 
are future-pointing unit timelike vectors and so the symmetry 
absent. Nevertheless there exist non-degenerate, 4-dimensional, stationary 
points as demonstrated in detail below. These correspond to 4-simplexes in 
Minkowski space in which all tetrahedral faces are spacelike. A new feature 
is that some of these faces are future-pointing and some are past pointing. 
However the scaling behaviour of these stationary points is exactly the same 
as for the SO (4) case. Similary the other generic cases of stationary points 
exist in the same way and exhibit exactly the same scaling as for SO (4). 
In particular, the one-dimensional stationary points again scale as q;~^ and 
again dominate the generic configuration in the asymptotic limit. 

A new feature in the case of non-generic spins is the possibility of station- 
ary points lying on the boundary of hyperbolic space. These will correspond 
to the possibihty of tetrahedral faces of the geometric 4-simplex which are 
null. Clearly null tetrahcdra cannot occur in the SO (4) case, where vanish- 
ing 3-volume implies that a tetrahedron is degenerate. However we do not 
analyse this interesting possibility in any further detail here. 

For most of the types of stationary point the equations are the same as 
for the SO (4) case, and so they are not repeated here. The one case which 
is analysed in detail is the non-degenerate case in the next section. 

4.3 Non-degenerate 4-dimensional configurations 

Consider configurations where the Xi are non-parallel and span 4 dimensions. 

Varying the action with respect to the y's we find that the action is 
stationary when yj^i is parallel or antiparallel to ^ki- Let — ±1 then 




(49) 



The action at the stationary points in the y's becomes 




(50) 



k<l 
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To examine the contribution from non-degenerate stationary points we need 
a Schlafli identity for Lorentzian 4-simplexes in which every 3-dimensional 
face is spacehke. The normal vectors to these faces are therefore timehke and 
we need a definition of the angle between them. 

Definition: Lorentzian angle jl3| The Lorentzian angle between two time- 
like unit vectors rik and rii, denoted by Qki, is defined in two distinct 

cases 

1. Interior Lorentzian angle: If one of and rii is future pointing 
and the is other past pointing then 9^^ > and is given by 

Qki = cosh'^ {nk.rii) (51) 



2. Exterior Lorentzian angle: If Uk and ni are both future or past 
pointing then G^/ < and is given by 

Qki = -cosh'^{-nk.ni) (52) 



Setting TUki = for exterior Lorentzian angles and rriki = 1 for interior 
Lorentzian angles we can combine these definitions into 

Qm = -i-ir'^' cosh-i(-(-l)™^'nfc.nO. (53) 

This implies 

nfc.n, = -(-1)™" cosMBh) 
for all k, I, as long as 9^^ is defined to be 0. 



Theorem 1 (Lorentzian Schlafli identity). For a Lorentzian ^-simplex 
with areas A^i and timelike normal vectors ui . . .n^ 

Y,AkidQu = ^ (54) 

k<l 

Proof: This follows the derivation of Schlafli identity in ^3] but using Lor- 
entzian instead of Euclidean geometry. Let cr be a 4-simplex with 
timelike normal vectors ui-.-Ur, and Lorentzian 4- volume \a\, af be 
the tetrahedron opposite vertex i with 3- volume Tj and afj the triangle 
of a opposite vertices i and j with area Aij. 

Stokes' theorem gives 

TkUk = 

k 
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and dotting this with ni imphes 

5^-(-ir- cosh(efc,)Tfc = 

k 

Differentiating, and contracting with the vector T; gives 

J2TkTismh\eki\deM = (55) 

kj^l 

However, one can show that 

4 

TfcT, sinh 16^;! = -Aki\cr\. (56) 
This follows from the equations 

\a\ = —Tkh 
II ^ k 

Ti = ^h'Aki 

where h and h' are the altitudes of the 4-simplex and tetrahedron T; 
respectively. These are related hj h = /I'sinh \ Qki\- 

Now equation simplifies to 

J2Akideki = 0. (57) 

k<l 

□ 

Consider a geometric 4-simplex with timelike normals ni...?7,5. Ro- 
tate the 4-simplex so that rii is at (1,0,0,0), n2 is in the plane spanned 
by (1, 0, 0, 0),(0, 1, 0, 0) and 77.3 is in the hyperplane spanned by (1,0,0,0), 
(0,1,0,0), (0,0,1,0). 

Now some of the normal vectors will be future pointing, and others will 
be past pointing. However the definition of the Xi vectors in equation ()47|1 
requires them to be future pointing. Take Xi = atrii where Oj = 1 if 
is future pointing and = — 1 if nj is past pointing. From the choice of 
alignment of the 4-simplex we always have Oi = 1. These definitions give 



(58) 
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Varying the action, and using a Lagrange multiplier A for the constraint 
dS" = ekiPkiauai d©M ^ X^Aki d©w (59) 

k<l kl 

Set A = ±1 to fix the scale of the 4 simplex. Then we can identify 

Pki = Ml (60) 

and 

efei = AafeC; (61) 

which fixes the e's up to an overall sign. 

The expected contribution from these configurations is an oscillating func- 

9 

tion multiplied by a~2. The stationary phase calculations are similar to the 
SO (4) case but now the sign of is determined by whether face k is future 
or past pointing. Not all = 1 or —1 is possible since the 4 simplex must 
have at least one face future pointing and past pointing. 
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5 Conclusions 

The stationary phase technique has been used to calculate asymptotic for- 
mulae for SO (4) relativistic spin networks. 

The asymptotic formula for the relativistic spin network tetrahedron 
agrees exactly with the square of the Ponzano-Regge formula whenever we 
have a tetrahedron which can be embedded into 3-dimensional Euclidean 
space. 

For the SO (4) relativistic spin network 4-simplex there are a number 
of types of stationary points in the asymptotic formula. For a generic set of 
spins the evaluation of the lOj-symbol is dominated in the asymptotic limit by 
degenerate one-dimensional geometries, whose contribution in general decays 
as a constant times o;"^. 

The stationary points for the SO (3, 1) relativistic spin network 4-simplex 
can be classified in the same way as for the SO (4) case. There exist non- 
degenerate, 4-dimensional, stationary points which correspond to 4-simplexes 
in Minkowski space in which all tetrahedral faces are spacelike. The sta- 
tionary points corresponding to degenerate one-dimensional geometries again 
scale as a~'^ and dominate the generic configuration in the asymptotic limit. 
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